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We give a criterion to test a non-biaxial behavior of infinitely thin hard platelets of 
D2h symmetry based upon the components of three order parameter tensors. We inves- 
tigated the nematic behavior of monodisperse infinitely thin rectangular hard platelet 
systems by using the criterion. Starting with a square platelet system, and we compared 
it with rectangular platelet systems of various aspect ratios. For each system, we per- 
formed equilibration runs by using isobaric Monte Carlo simulations. Each system did 
not show a biaxial nematic behavior but a uniaxial nematic one, despite of the shape 
anisotropy of those platelets. The relationship between effective diameters by simula- 
tions and theoretical effective diameters of the above systems was also determined. 

KEYWORDS: isotropic-nematic phase transition, isobaric Monte Carlo simulations, coexistence number 
density, effective diameters 

1. Introduction 

Liquid crystals are molecules with some preferred orders without any positional or- 
der. Uniaxial nematic liquid crystals have one orientational order vector, referred to 
as a director n. When nematic liquid crystals have not only the first director n but 
also a secondary director m which is orthogonal to n (FIG. 1.), the nematic liquid 
crystals are said to be biaxial. They are equipped with a 10-100 times quicker response 
time under an external field than uniaxial ones due to their anisotropic shape. 1 (See 
FIG. 2.) In 1970, based on the molecular field theoretical calculations, Freiser pre- 
dicted that anisotropic-shaped molecules should undergo a biaxial nematic phase. 2 The 
first lyotropic (concentration dependent) biaxial nematic liquid crystals were synthe- 
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sized by Yu and Saupe in 1980. 3 Twenty-four years later, the discovery of a system of 
thermotropic (temperature dependent) biaxial nematic liquid crystals, based on an oxa- 
diazole unit, was announced by Maden, Dingemans et. al., 4 ' 5 but later, it was realized 
that, in fact, the system consisted of a uniaxial plane of smectic C liquid crystals. 6 Since 
then, new experimental methods have been developed to detect biaxial liquid crystals, 
such as 2 H (Deuterium) nuclear magnetic resonance 7 and IR (Infrared) 8 spectroscopies, 
but it remains a challenge to study the biaxial properties of nematic liquid crystals in 
the laboratory. 6 ' 9 ' 10 




X 



Fig. 1. The first and secondary directors n and m by which the third 1 = n x m is defined auto- 
matically, n, m, 1 are called laboratory coordinate axes. On the other hand, molecular fixed coordinate 
axes are indicated by x, y, and z. 




Fig. 2. A model of anisotropic shaped liquid crystal molecules. The motion b is 10-100 times quicker 
than the motion a. 9 

The first theoretical study of phase transitions in liquid crystals were carried out by 
Onsager in 1949 using hard rod and disc models. 11 This work was extended with the 
advent of a Monte Carlo simulation technique 12 by Eppenga and Frenkel in 1984 who 
modeled the isotropic-nematic phase transition of the system of monodisperse infinitely 
thin circular platelets by excluded volume effects. 13 Bates and Frenkel detected the first 
order nematic-columnar phase in the system of the circular platelets 14 with a scaling ar- 
gument introduced by Bolhuis and Frenkel. 15 Bates studied the isotropic-nematic phase 
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transition of monodisperse infinitely thin circular and non-circular (regular hexagon, 
pentagon, square, triangle, ellipse and rectangular of a width to height aspect ratio of 
1:2) platelet systems by using grand canonical Monte Carlo simulations. He pointed 
out that the shape of a platelet is important to determine its isotropic-nematic phase 
transition. 16 

The plan of this paper is as follows. First, we present a set of sufficient conditions that 
one can use to test non-biaxiality (See FIG. 3) in a system of infinitely thin hard platelet- 
like molecules of D 2 h symmetry. Next, as a concrete example, we study the behavior 
of square and various rectangular platelet systems which have not been fully examined 
in the literature, by isobaric Monte Carlo computer simulations. Using the sufficient 
criterion, we check if those systems are non-biaxial. The third, we investigate how the 
phase transition pressure of each platelet system depends on the shape anisotropy of 
platelets, starting from square platelets to elongated rectangular ones. 




n, m are both determined. 

v J 



Fig. 3. A flow chart for uniaxial, non-biaxial and biaxial tests. The non-biaxial test is shown by a 
shaded area. 

2. Criterion for non-biaxiality 

For liquid crystal molecules of arbitrary shape, a generalized order parameter tensor 
S is a real symmetric (and therefore diagonalizable) matrix quantity which measures a 
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phase transition. Its components are given by 

S% = ^.(^r^-^a^, (1) 

= ^(a* • r »)( a j ■ **) ~ \ s ij s u^, (2) 

where (•) is an ensemble average, 5 is the Kronecker delta function, a^j are unit vectors 
along three laboratory coordinate axes 1, m, n, and r^ u refers to the unit vector parallel 
to three molecular fixed coordinate axes x, y,z (FIG. 1), and I is the 3x3 identity 
matrix. 17 All i, j, /i, v take the values 1, 2, and 3. For the components of a nematic order 
parameter tensor, we set ji = u, and then we obtain: 

^ = (^(a l -r,)(a r r M )-^,y (3) 
In the system of iV molecules, Eq. (3) is computed as 

S % = ( a7 E {l^ ■ • r» - (4) 

Here, the index k attached to r M runs over all molecules in the system. 

The largest absolute value of eigenvalues of the diagonalized matrix is called a 
uniaxial order parameter S£. n , and its corresponding eigenvector is the first director n. 
We simplify the notation by writing as S zz = S z n z n , S yy = S"" and S xx = S" x • We can 
also write them in terms of Saupe parameters, S = S zz , D = S'xx — S'yy. 18 To study a 
uniaxial behavior of anisotropic molecules, it has been often the case that one evaluates 
D by taking the difference of S xx and S'yy; however, we keep the raw values of S xx and 
S'yy as they are. It is because they can tell us a non-biaxial nematic behavior of infinitely 
thin anisotropic molecules of D 2 h symmetry. More precisely speaking, assuming that 
those anisotropic molecules show a uniaxial nematic phase (S zz > 0.7), if we have S yy ~ 
S xx ~ —0.5, then they are in a non-biaxial phase. Supposing in stead that they were in a 
biaxial phase, we would have S zz > 0.7, S yy > 0.7, S xx > 0.7. Without any assumptions 
such as infinitely thin and D 2 h symmetry, we need to confirm S zz > 0.7, C > to say 
that arbitrary anisotropic molecules are in a biaxial nematic phase. Here, C is another 
Saupe parameter, and it is defined by C = S" x + Sy^ m — S yy — S^ m . 18 All numerical 
values used here are estimated values to describe a typical nematic phase. 

The magnitude of the difference between the absolute values of the second and third 
largest eigenvalues S ZZj 2, S ZZj 3 of the matrix S z {, which we shall call a rocking parameter, 
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is denoted by R zz = ||£ zz ,2| — |<5 Z z,3||- This measures to what degree a system shows 
a rocking-chair like nematic behavior with two orthogonal preferred orientations. In 
general, the eigenvectors corresponding to the second and third largest eigenvalues are 
not always the second and the third director m and 1 (= n x m). The other rocking 
parameters with respect to y , x axies R yy = \ \ S yy ^ \ — \ S yy ,3 \ | and i? xx = 1 1 S'xx^ | — | 5xx,3 1 1 
are defined similarly. Rocking parameters are relatively small; nevertheless, if a system 
is in a non-biaxial nematic phase, we have the following relation among them: R zz < 
{i? yy , i? xx }. Dingemans et. al. called R zz as a biaxial molecular order in contrast to a 
phase biaxiality rj = S™ m — 5^. 5 To the best of our knowledge, the meaning of R yy 
and i? xx has not been investigated yet. 

In summary, we give the above non-biaxial criteria in a comparison with uniaxial 
and biaxial conditions in Table I. 

Table I. The list of phases and the set of parameters to be determined. We cannot determine that 
a system is non-biaxial only by checking the condition for being a uniaxial phase. The information of 
5 V v an d <Sxx are needed. 



Phase 


Parameter 


Condition 


Uniaxial 






<Szz 


> 0.7 


Non-biaxial 


SzZ 7 




<Szz 


> 0.7, Syy W Sx X W -0.5 


Biaxial 


S zx , 


c 


*Szz 


> 0.7, C > 



3. Models 

Monodisperse infinitely thin rectangular (of aspect ratios 1:1, 1:2, 1:2.5, 1:3, 1:3.5, 
1:4, 1:5) hard platelet systems (FIG. 4) are used. Previously, we examined the isotropic- 
nematic phase transition of an N = 120 square (1:1) platelet system. 19 This time, each 
system is composed of iV = 480 platelets of equal area A. We set the areas of square 
and rectangular platelets to be the same. The reduced pressure is set by p* = D 3 /3p, 
where the unit length D is defined by D = (v^4)/2, f3 is the inverse temperature, and p 
is a pressure. We performed equilibration runs using isobaric Monte Carlo simulations. 
The initial configuration of each system was isotropic. We isotropically compressed it 
at reduced (dimensionless) pressures p* ranged in p* = 0.1-2.0. An excluded volume 
repulsion drives the phase transition of each system. We compute a number density 
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p = Na 3 /V, where N is the number of platelets, a = ^/AA/tt is a scaled diameter, and 
V is the volume of a periodic simulation box. About up to 50000 Monte Carlo steps are 
typically used until equilibrium is reached. Thermodynamic quantities were calculated 
with up to 10000 Monte Carlo steps. 



Fig. 4. Square (1:1) and rectangular platelets of equal area with aspect ratios 1:2 and 1:5 are shown 
as examples. The z-axis is coming out of the paper where the x- and y-axes cross. 



4. Numerical results and Discussion 

The plots of reduced pressure p* vs. the number densities p and S xx vs. p from 
an isotropic phase to a nematic phase are shown in FIG. 5. We computed a phase 
transition point by taking the center of the steepest joined segment of two data points. 
In our previous study, the coexistence number density p t (the subscript t means a 
transition number density) for 120 square platelets was determined by using the phase 
space multi-histogram method (PSMH), and it was p t = 3.87 ± 0.23 at p* = 1.13. 19 
In contrast, taking the center of the steepest joined linear segment of two data points, 
we found that the coexistence number density of the 120 square platelet system was 
p t = 3.60 ± 0.11 at p* = 1.13. The results from both methods corresponded well within 
errors. As an aspect ratio increased, a phase transition pressure (or coexistence number 
density) decreased. In Table II and FIG. 6, we summarized phase transition reduced 
pressures and coexistence densities and the phase diagram for platelets of various aspect 
ratios. 

Before going into the details of our further results, we shall recall the definition of 
effective diameters. Let us denote the number densities of the system of rectangular 
platelets by p ai . The lower index i adjacent to a indicates the aspect ratio of a rectan- 
gular platelet a written in the form i = i/1 instead of 1 : i. As Bates defined an effective 
diameter to be the cube root of the ratio of the isotropic-nematic coexistence number 
density for a non-circular platelet system to a circular platelet system, 16 we analogously 
define by a cS = -sj p a i/p a i- Here, p a x, p ai are coexistence densities, and strictly written, 



6/16 



J. Phys. Soc. Jpn. 



FULL PAPERS 





Fig. 5. (a) Pressure dependance of the number densities p. p* means a reduced pressure, (b) 5 ZZ vs. 
p. We started from isotropic and compressed each system to nematic. Throughout this paper, unless 
otherwise stated, the segment lines drawn between data points are intended to guide the eye, and error 
bars are the standard deviation of the mean. 
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Table II. Phase transition reduced pressures (p*) and coexistence number densities (p t ) of platelets 
of various aspect ratios. 



Aspect ratio 


P* 


Pt 




1. lO 


O fif) J. All 


1:2 


l.n 


3.52 ± 0.07 


1:2.5 


1.05 


3.34 ± 0.07 


1:3 


1.03 


3.27 ± 0.05 


1:3.5 


0.98 


3.11 ± 0.06 


1:4 


0.95 


2.97 ± 0.04 


1:5 


0.88 


2.83 ± 0.06 



3.8 



3.6 



3.4 



Cl 3.2 



3.0 



2.8 



2.6 
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Fig. 6. Phase diagram for infinitely thin rectangular hard platelets. The aspect ratio i of the platelets 
is in the form i = i/1 instead of 1 : i. The dotted line drawn through the data is the phase boundary 
of isotropic and non-biaxial nematic. 



it might be denoted as pt, a %- For simplicity, we adapt the notation with omitting the 
subscript t indicating a phase transition as before. Following Kihara's work, 20 Bates 
computed the orientation ensemble average excluded volume (denoted by (V*f ) in this 
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Fig. 7. Effective diameters from simulations vs. theoretically calculated effective diameters. The 
solid line indicates a c ff — <r cS , and the obtained data was fitted by the dotted line 1.07<7 e ff — 0.06 = cr off . 



paper) of a pair of infinitely thin platelets of arbitrary convex shape and found that 
(y*f) = \AP ai , where A is the area of a platelet ai, and P ai is its perimeter. 16 In the 
same spirit as Bates', taking the cube root of the ratio of excluded volumes, we de- 
fine a theoretical effective diameter a' cS = \/ (V^f ) / (V^f ). It is of use to check whether 
°cff — °cfr holds or not since if so, one can predict systems of other aspect ratios by com- 
puting the orientation average excluded volume of a pair of platelets without performing 
further simulations. In an earlier study, it was expected that a cS = a' cS will not hold for 
rectangular platelets as their aspect ratio is bigger than 1:2 since circular platelets were 
taken instead of square platelets (See FIG. 4.) 16 We see clearly in FIG. 7 that a cS = a' cS 
did not hold, but there observed that the measured effective diameters by simulations 
were smaller than the theoretical ones. The fitted line had its slope a' eG /a eS ~ 1.07. 

FIG. 8 displays p* vs. the order parameters S zz , S yy , S xx and the rocking parameter 
with respect to z-axis R zz of the systems of rectangular platelets of aspect ratios 1:1, 1:3 
and 1:5. (Those for other aspect ratios are not shown here.) Following the S zz plot, the 
first order phase transition pressure gradually became low as the aspect ratio increased. 
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The non-biaxial (uniaxial) nematic phase for square and rectangular platelet systems 
was detected from S zz > 0.7, S yy « S xx « —0.5. We found a low rocking parameter 
R zz < 0.25 for every shape. A physical meaning of the values of R yy (or R xx ) is that 
platelets swing like a rocking chair as they rotate around each z-axis. In FIG. 9, we show 
p* vs. R zz , R yy , R xx of the systems, and we found the relationship R zz < {R yy , R xx }- We 
also found that S zz > 0.7 and S yy < (or S xx < 0) from FIG. 10. Among the systems, 
fluctuations (variations in S zz ) of lower aspect ratios were bigger than higher ones near 
their phase transition points. (See FIG. 11) These indicated that the more elongated 
platelets showed a uniaxial nematic phase as they were rotating around the z-axis in the 
more stable manner, and in addition, they swung along the y- and x-axes. (See FIG. 12.) 
This could be observed simply by computing the components of order parameter tensors 
with respect to the z, y, and x-axes. If it were in a biaxial nematic phase, we would 
expect to have S zz > 0.7, S yy > 0.7, S xx > 0.7, and R zz « R xx « R xx 0. 

5. Conclusions 

We introduced a non-biaxiality criterion for infinitely thin anisotropic molecules 
of D 2 h symmetry. This is a test for uniaxiality, but it cannot be applied to prove 
biaxiality. To confirm biaxiality, we would need to add additional constraints on the 
C parameter which can be determined 18 with additional simulations. Our criteria for 
non-biaxiality offers an efficient means to rule out biaxiality. We studied the isotropic- 
nematic behavior of monodisperse infinitely thin square and rectangular hard platelet 
model systems by isobaric Monte Carlo simulations. Although one might expect that 
the shape anisotropy would generate a biaxial nematic phase, we demonstrated the 
shape anisotropy of the platelet systems, from square platelets to rectangular ones, did 
not induce a biaxial nematic behavior by calculating components of order parameter 
tensors. Those platelets were in a non-biaxial (uniaxial) nematic phase at a high reduced 
pressure. We observed that they were rotating around the z-axis as they swung along 
the y- and x-axes. The sufficient criterion to check systems of monodisperse infinitely 
thin liquid crystal molecules of D 2 h symmetry being non-biaxial was presented. The 
relationship seen in effective diameters of those systems was determined. 
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Fig. 8. Pressure dependance of the order parameters 5 ZZ , S yy , S xx and the rocking parameter R z . 

in the systems of 480 rectangular platelets of aspect ratios (top) 1:1, (middle) 1:3 and (bottom) 1:5. 
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Fig. 10. S zz vs. S yy or 5 XX . (a) 1:1, (b) 1:3 and (c) 1:5. 
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Fig. 11. (a) Variances (abbreviated by Var.) of 5 ZZ , written more precisely (5S% Z ), with p* /pi for 
various aspect ratios. Here, p* t is the reduced pressure of the isotropic-nematic phase transition, (b) 
The magnified plot of (a). 
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